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Effects of wall cooling on the receptivity process induced by the interaction of slow 
acoustic disturbances in the tree-stream are numerically investigated for a boundary layer 
flow over a 5-degrees straight cone. The tree-stream Mach number is 6.0 and the Reynolds 
number is 7.8x10® /ft. Both the steady and unsteady solutions are obtained by solving the full 
N avie r-Stokes equations using 5 th -order accurate weighted essentially non-oscillatory 
(WENO) scheme for space discretization and using 3 rd -order total variation diminishing 
(TVD) Runge-Kutta scheme for time integration. Computations are performed for a cone 
with nose radius of 0.001 inch for adiabatic wall temperature (T aw ), 0.75*T aw , 0.5*T aw , 
0.40*T aw , 0.30*T aw , and 0.20*T aw . Once the mean flow field is computed, disturbances are 
introduced at the upstream end of the computational domain. Generation of instability 
waves from leading edge region and receptivity of boundary layer to slow acoustic waves are 
investigated. Computations showed that wall cooling has strong stabilization effect on the 
first mode disturbances as was observed in the experiments. Transition location moved to 
upstream when wall cooling was applied It is also found that the boundary layer is much 
more receptive to fast acoustic wave (by almost a factor of 50). When simulations performed 
using the same forcing frequency growth of the second mode disturbances are delayed with 
wall cooling and they attained values two times higher than that of adiabatic case. In 
0.20*T aw case the transition Reynolds number is doubled compared to adiabatic conditions. 
The receptivity coefficient for adiabatic wall case (804 °R) is 1.5225 and for highly cooled 
cones (241, and 161 °R); they are in the order of 10" 3 . 


I. Introduction 

T ransition process from laminar to turbulent flow is still an important challenge even after years of research. 

Accurate prediction of transition location is vital for the design of hypersonic vehicles. Because transition 
controls important quantities such as aerodynamic drag, heat transfer and other boundary layer parameters. In 
hypersonic boundary layers one important parameter is the wall temperature. Wall cooling would be expected to 
stabilize first mode disturbances while destabilizing the second mode. The effects of cooling on transition have been 
studied experimentally and numerically by many researchers. 

Lees 1 predicted that cooling the surface would stabilize the boundary layer. Later, Mack’s 2 results showed that 
the first mode was stabilized by cooling; however the higher modes were destabilized by this process. Experiments 
of Demetriades 3 and Lysenko 4 confirmed that cooling the wall increased the growth rates of the second mode 


Graduate Research Assistant, Student Member AIAA 
: Research Scientist, Member AIAA 
* Professor & Eminent Scholar, Associate Fellow AIAA 


1 

American Institute of Aeronautics and Astronautics 



disturbances and reduced the transition Reynolds number. Stetson 5 investigated the effects of surface temperature on 
the stability of the laminar boundary layer experimentally and the results also verified linear stability theory. 
Balakumar and Malik 6 computed the parametric effects of the pressure gradient and wall cooling on the stability of 
the flow over cone using the quasi-parallel, compressible linear stability e Mallk code. 

Receptivity 7 is a process by which free-stream or wall-induced disturbances enter the laminar boundary layer and 
generate unstable waves 8 . The receptivity process is poorly understood at hypersonic speeds. The understanding of 
this phenomenon is of great importance because receptivity connects the amplitude of the free-stream disturbances 
and initial amplitude of the unstable waves 9 . Recent experimental 9 " 13 , theoretical 14 ' 19 and computational 20 ' 35 studies 
increased our understanding about receptivity mechanism. However it still remains as a challenging problem with 
practical importance. 

The objectives of this work are to estimate the destabilizing effect of the wall cooling on the hypersonic 
boundary layers over a blunt cone and to calculate the receptivity coefficient of the instability waves generated near 
the leading edge. To do so simulations are performed at different wall temperatures T aw (adiabatic wall temperature), 
0.75*T aw , 0.50*T aw , 0.40*T aw , 0.30*T aw and 0.20*T aw for leading edge nose radius R n =0.001 inch and at a unit 
Reynolds number of 7.8*1 0 6 /ft for a 5-degree half angle cone. 
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Applying same procedure to shear stresses and heat fluxes we get new term without singularity problem at y=0. 
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Here ( X f y) are the two-dimensional and (X, y, 6) are the axisymmetric coordinates and ( U , V) are the 
corresponding velocity components, p is the density, p is the pressure. E is the total energy given by 
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Here e is the molecular internal energy and T is the temperature. 

The viscosity (//) is computed using Sutherland’s law and the coefficient of conductivity (A) is given in terms of 
Prandtl number {Pi). The variables p , p, T and velocity are non-dimensionalized by their corresponding reference 

variables p^ /?«,, 7" OT and V RT respectively. The reference value for length is computed by -\JvX Q / U :f , where Xq is 

the location of the beginning of the computational domain in the stream wise direction. For the computation, the 
equations are transformed from physical coordinate system {X, ]/) to the computational curvilinear coordinate system 
(E,,r|) in a conservative manner and the governing equations become 
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I II. Solution Algorithm 

The governing equations are solved using 5 th order accurate weighted essentially non-oscillatory (WENO) 
scheme for space discretization and using 3 rd order total-variation-diminishing (TVD) Runge-Kutta scheme for time 
integration. These methods are suitable in flows with discontinuities or high gradient regions. These schemes solve 
the governing equations discretely in a uniform structured computational domain in which flow properties are 
known at the grid nodes. WENO scheme approximate the spatial derivatives in a given direction to a higher order at 
the nodes, using the neighboring nodal values in that direction. TVD-RK scheme integrates the resulting equations 
in time to get the point values as a function of time. Since the spatial derivatives are independent of the coordinate 
directions, the method can easily add other dimensions. It is well known that approximating a discontinuous 
function by a higher order (two or more) polynomial generally introduces oscillatory behavior near the 
discontinuity, and this oscillation increases with the order of the approximation. The essentially non oscillatory 
(ENO) and the improvement of these WENO methods are developed to keep the higher order approximations in the 
smooth regions and to eliminate or suppress the oscillatory behavior near the discontinuities. They are achieved by 
systematically adopting or selecting the stencils based on the smoothness of the function, which is being 
approximated. Shu 36 explains the WENO and the TVD methods and the formulas. Atkins 37 gives the application of 
ENO method to the N-S equations. Balakumar 38 describes in detail the solution method implemented in this 
computation. 

Four different types of boundary conditions are used in steady and unsteady calculations. At the outflow 
boundary, extrapolation boundary condition is used. At the wall, the simulation uses viscous conditions for the 
velocities and adiabatic or isothermal wall conditions, and it computes density from the continuity equation. In the 
mean flow computations, the simulation prescribes the free-stream values at the outer boundary, which lies outside 
the bow shock. The blunt cone is assumed to align with the free stream. In the unsteady computations, it 
superimposes the acoustic perturbations to the uniform mean flow at the upper boundary. The procedure is to first 
compute the steady mean flow by performing unsteady computations using a variable time step until the maximum 
residual reaches a small value ~10" n . These computations use a CFL number of 0.2 for adiabatic case and 0.1 for 
isothermal cases. The next step is to introduce unsteady disturbances at the upper boundary of the computational 
domain and to perform time accurate computations to investigate the interaction and evolution of these disturbances 
to downstream. Since we use very fine spatial grid to resolve the leading edge region, these computations require 
very small time step which is taken as the minimum time step allowable for aforementioned CFL numbers. 
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The grid is generated using analytical formulae. The grid stretches in the rj direction close to the wall and is 
uniform outside of the boundary layer. In the E, direction, the grid is symmetric about the leading edge and very fine 
near the nose and is uniform in the flat region. The outer boundary that lies outside the shock follows a parabola so 
that the boundary layer growth could be captured accurately. The computational domain extends from x=-0.015 to 
20.0 inches in the axial direction in this computation depending on the wall temperature. Calculations were 
performed using a grid which has 32 blocks and each block has 127x425 grid points. Computational domain has 
approximately 2 million grid points. Due to the very fine grid requirement near the nose, the allowable time step is 
very small and the computations become very expensive to simulate the unsteady computations in the entire domain. 

The acoustic field that impinges on the outer boundary is taken to be in the following form. 


fj = Real ^pd aacX± ' eacy ~ lmt j 

Here a 3C , s ac are the a; y wave numbers, respectively, of the acoustic wave and 0 ) is the corresponding frequency 
of the acoustic disturbance. The incident angle Oof the acoustic wave is defined as 

0 = tan -1 

a ac 

and in this study computations are performed for zero incidence angles. 


IV. Results and Discussion 

Computations are performed for hypersonic flow at a free stream Mach number of 6.0 over a 5 -degree half-angle 
cone with blunt leading edge, R n =0.001 in, for different wall temperatures to investigate the effects of wall cooling 
on hypersonic boundary layer receptivity due to acoustic disturbances in slow mode. The flow parameters are given 
in Table 1 and the boundary layer edge conditions for a sharp cone are given in Table 2. Figure 1 shows the 
schematic diagram of the computational set up. The nose region of the cone is modeled as a circle. Simulations are 
performed for wall temperatures T w = Adiabatic wall (T aw ), 0.75*T aw , 0.50*T aw , 0.40*T aw , 0.30*T aw , and 0.20*T aw . 
Different cases are summarized in Table 3. 


Table 1. Flow parameters for the wind tunnel model (Horvath et. al. 39 ) 


Free stream 

Mach number 
Reynolds number 
Density 
Velocity 

Reservoir Pressure 
Reservoir Temperature 
Wall temperature 
Prandtl number 
Ratio of specific heats 
Length scale 


M ro = 6.0 
Re ro = 7.8xl0 6 / ft 

Poo = 7.059xl0' 3 lbm/ft 3 

Uoo = 3140.21 ft/s 

P 0 = 475 psi 

T 0 = 475 °F 

Adiabatic and Isothermal conditions 
Pr = 0.72 

y = 1.4 

^v m x 0 / U x ~ 7.308 x10 3 ft 


Table 2. Parameters along the cone surface. 


(Var.) edae / (Var.)oo 

Sharp Cone 

Mach Number 

0.932 

Pressure Ratio 

1.560 

Density Ratio 

1.372 

Temperature Ratio 

1.137* 

* Adiabatic case. 
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Table 3. Wall temperatures and bow shock standing distances. 


T wall I 

' Ad.wall 

r ^ wall / 

C/ 0 - j 

1 00 

T W all(° R) 

8 (x 10 4 in.) 

1.00 

7.052 

804 

2.50 

0.75 

5.289 

603 

2.39 

0.50 

3.526 

402 

2.28 

0.40 

2.8208 

322 

2.21 

0.30 

2.1156 

241 

2.15 

0.20 

1.4104 

161 

2.13 


T m = 113.98 °R. 


Linear Stability 

Linear stability computations done for flow over an axisymmetric cone at a free stream Mach number of M=6 at 
different wall temperature conditions Tw/T aw =1.0, 0.50 and 0.20. Figures 2(a), (b) and (c) depict the results in {Re, 
F ), {Re, a) and ( Re,C r ) planes for two-dimensional disturbances respectively. It is seen that with the wall cooling (1) 
the first mode is stabilized (2) the unstable frequencies increase (3) the unstable wave number increases and (4) the 
range of the phase speeds of the unstable waves narrows. 

Figure 3 (a) -(d) show the N-Factor curves for decreasing frequencies at different wall temperature conditions. 
Most amplified frequencies for wall temperatures TJT^I.O, 0.75, 0.50 and 0.20 are found as 0.8x1 (J 4 , 0.9x1 (J 4 , 
1.1x1 Of 4 , and 1.7x1 Of 4 for N=10. Here the variables are non-dimensionalized by the variables at the edge of 
boundary layer. To obtain the variables non-dimensionalized by the free stream values as given in Table 1, the 
variables in this section should be multiplied by the appropriate factors from Table 2. The frequency variable ^has 
to be multiplied by 1. 174 to obtain the values in terms of free stream values. 

The neutral stability curve clearly shows the unstable first and the second mode regions for adiabatic wall 
temperature but when cooling applied the unstable first mode region disappears. The first mode and the second 
mode neutral stability curves merge at a Reynolds number of Re= 1600 fox the cone. 

Mean Flow 

Validation studies of axisymmetric code presented in a previous study 33 . Figure 4 shows the mean flow 
temperature contours computed using the WENO code. The figures 5(a)-(d) show the results for the 5^degree half- 
angle cone at different wall temperatures T w =AdiabatiC wall (T aw ), 0.75* T aw , 0.50 *T am and 0.20*T aw . In the 
adiabatic wall case bow shock generated a high temperature region. This region convected to downstream over cone 
wall. For cooled wall cases the high temperature region is trapped between the bow shock and nose part of the cone. 
Bow shock standing distance reduced %15 when wall temperature reduced to 0.20*T aw (Table 3). Figure 5(a) and 
5(b) depict the density contours over adiabatic cone and cooled cone ( T w =0.20*T aw ) cases. Maximum non- 
dimensional density occurred on the stagnation point for adiabatic and cooled conditions are 6 and 33 respectively. 
Over the flat part of the cone density contours looks similar. High density region is increased when cooling applied 
but it did not extend to flat part. Figure 6 compares the Mach contours of the aforementioned wall conditions. Wall 
cooling decreased the oblique shock angle and compressed the Mach layers to wall. Figure 7(a) shows the wall to 
free stream temperature ratio, Go , along the cone surface for different wall cooling conditions and Fig. 7(b) shows 
the variation of boundary layer edge Mach number, M e , along cone surface. Edge Mach number and temperature for 
adiabatic case at x= 13.96 in. are M e =5.575 and T e = 129.53 °R and for 80% cooled case at x= 19.58 in. are M e =5.581 
and T e = 129.26 °R. 

Density profiles at different x locations are plotted in Figures 8(a)-(d) in similarity coordinates for adiabatic 
wall and cooled walls {0.50 *T aw , 0.30 *T aw , and 0.20 *T av )). At X=0.2 in. density profiles converges to a profile which 
stays the same until the end of the computational domain. Wall cooling decreased boundary layer thickness from 
rj=14.6 ( T aw ) to rj=8.19 ( 0.20*T aw ) and also increased the density on wall from 0.16 ( T aw ) to 0.80 ( 0.20*T aw ). 
Density increase on the wall changed the characteristic of density profiles. In adiabatic case minimum density 
occurred on cone wall and steadily increased away from it. However, for cooled wall case ( T w =0.20*T ay[ ) density on 
the wall (p r}= o =0.80) first started to decrease away from the wall until rj=2 (p r?=2 =0.42). After making a minimum 
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at this point density changed its character and increased as in the adiabatic wall case. Even this density profiles 
seems to have a low density region their minimum density at rj=2 (p v= 2 =0.42) is 2.6 times larger than the minimum 
of adiabatic case (p v=0 =0.16). As a result wall cooling reduced the boundary layer thickness and increased the 
density inside it. 

Figures 9(a)-(f) show the temperature profiles at different axial locations for different cooling cases. The 
compressible Blasius similarity profiles are also included for comparison. In adiabatic wall case at X=0.2 in. 
calculated temperature profiles perfectly matched with similarity solution (Fig. 9(a)). As was observed in the 
density graphs, temperature profiles are also converged to the same profile at X=0.2 in. and stayed same until the end 
of the computational domain. Wall cooling decreased thermal boundary layer thickness also. For highest cooling 
case ( T w =0.20*T aw ) wall temperature was 161 °R and it increased away from wall until rj=2 (T v=2 =305 °R). After 
making a maximum at this point temperature of the flow decreased as in the adiabatic case to 1 29. 26 °R. The 
difference between the similarity profiles and simulation results comes from the nose bluntness, R n = 0.001 in. Figure 
1 0(a) and (b) show the mean flow density profiles at different axial locations in physical coordinates for adiabatic 
( 804 °R) and highest wall cooling case ( 161 °R) respectively. Minimum density observed at y =0.01 in. normal to 
wall. Figure 11(a) and (b) depicts the mean flow temperature profiles at different axial locations in physical 
coordinates for adiabatic ( 804 °R) and cooled wall ( 161 °R) conditions. 

Interaction of Slow Acoustic Waves with the Boundary Layer 

After the mean flow is computed two dimensional slow acoustic disturbances are introduced at the outer 
computational boundary and the time accurate simulations are performed. Unsteady simulation results are presented 
for the cases T w = T aw , 0.75 *T aw , 0.50*T aw , and 0.20*T aw (Slow and Fast waves) at the most amplified frequencies 
0.96x1 Of 4 , 1.08x1 (I 4 , 1.32x1 Of 4 , and 2.0x1 (J 4 . These frequencies correspond to 375, 421, 515, and 780 kHz 
respectively. These frequencies give maximum amplification within the computational domain. To remain in the 
linear regime, the amplitude of the forcing freestream acoustic waves is given a small value of 2x1 Cf 5 for the first 
three cases and 2x1 (J 6 for 0. 20*T aw cdL$Q. Even with these small amplitudes, nonlinearity starts to develop near the 
end of computational domain for adiabatic case. 

Figures 12(a)-(e) show the evolution of the wall pressure fluctuations for aforementioned cases in linear scale 
while Figures 13(a)-(e) depict the results in a log scale. The figures clearly show the generation and the eventual 
exponential growth of the instability waves inside the boundary layer for adiabatic and cooled wall conditions. For 
the adiabatic wall conditions, Fig. 12(a) and Fig. 13(a), the disturbances grow from the leading edge and reach large 
amplitude levels of 0. 50 near the predicted transition onset point at X=19in. The slow wave whose wavelength is 
closer to the wavelength of the instability wave transforms into instability wave smoothly near the nose region. The 
parallel linear computations show that the first mode amplifies starting from leading edge weakly up to X ~ 15 inch 
for this frequency. 

Figure 13(b) shows the pressure fluctuations for T w =0.75*T aw m log scale. It is evident that wall cooling reduced 
the growth of the first mode until X=6 in. then growth of the second mode is observed. Because of wall cooling 
maximum amplitude only reached to 0.734 which is in the same order of adiabatic case. When more wall cooling 
( T w =0.50 T av J) is applied stabilization of the first mode become obvious. In Figure 13(c) the neutral point appeared 
at X=5.55 in. with an amplitude of 6. 18x1 O' 6 . Maximum amplitude for this case is 0.041 at X= 10.73 in. The neutral 
points are not discernable for adiabatic and T w = 0. 75 *T aw but it appeared for T w = 0. 50 *T aw and stayed in the picture 
for cooler cases. 

Figure 13(d) and (e) show the wall pressure fluctuations for slow and fast waves respectively at T w =0.20*T aw . In 
these figures the neutral points are located at 4.94 in. (Slow Wave) and 5.62 in. (Fast Wave) with the amplitudes of 
2x1 (J 7 8 and /^/^respectively. Maximum wall pressure amplitudes 0.00023 for slow wave and 0.0029 for fast wave 
are observed at x=7.6 in for both cases. Because of the difference in the initial amplitude of the acoustic disturbances 
the maximum amplitudes for adiabatic wall, 0.75 *7"^ and 0.50 *7"^ need to be divided by 10 and they become 0.05, 
0.0734, and 0.0041 respectively. While the maximum amplitude of 0.20 *T aw (Slow Wave) simulation, 0.00023, is 
almost 20 times less than 0.50 *T aw simulation, maximum amplitude of 0.20 *T aw (Fast Wave) simulation is on the 
same order. 

Also it is interesting to observe that the first mode is not stabilized for fast wave case in Figure 13(e) while 
Figures 13(a)-(d) clearly show the stabilization effect of wall cooling on the first mode for slow wave cases. 
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Simulated transition locations for T w =T aw , 0.75 *T aw , 0.50*T aw , and 0.20*T aw (Slow and Fast waves) wall 
temperature conditions are 19in, 13 in, 10.7 in, and 7.6 in. respectively. From these results one can conclude that wall 
cooling is destabilizing the boundary layer and transition points are moving to upstream. 

Another set of simulations run for T w = T aw , 0.75 *T am 0.50*T aw , 0.40*T aw , 0.30*T aw , and 0.20*T aw using the same 
forcing disturbance frequency, F= 1.2x1 (7 4 , in slow wave mode to observe the effects of wall cooling on 
development of instability waves in the same disturbance environment. Figure 14 and 15 shows the unsteady wall 
pressure fluctuations for aforementioned simulations. 

Figure 15(b) shows the pressure fluctuations for T w =0.50*T aw in log scale. It is evident that wall cooling reduced 
the growth of first mode until X=6.96 in where the amplitude is 6.5x1 5 6 then growth of second mode observed. 
Because of the wall cooling maximum amplitude only reached to 0.0035 which is 35 times less than of the adiabatic 
case. Stabilization of the first mode become apparent with wall cooling as was seen in neutral stability diagrams 
(Figure 2(a)). Figure 15(c) shows the exponential decrease of amplitude of wall pressure fluctuations for 
T w =0.40*T aw simulation until the neutral point located at X=8.46 in. with an amplitude of 9.95x1 5 7 . Maximum 
amplitude for this case is 0. 129 occurred at X= 14.62 in. More cooling increased the maximum amplitude level to the 
same value of adiabatic case. It clearly shows the destabilization effect of wall cooling on second mode. Finally for 
T w =0.20*T aw simulation the maximum amplitude increased more than two times and reached to 0.301 at the end of 
the conroutational domain X=19.9in. Further cooling moved the neutral point to x=13.4 in where the amplitude is 
4.7x1 5 8 . In these simulations neutral point was not discernable for adiabatic wall and T w = 0.75 *T aw but it appeared 
for T w = 0. 50 *T aw and cooler wall temperatures because of the stabilization of the first mode (Slow Wave) with wall 
cooling. Transition locations computed for T w =T aw , 0.50*T aw , 0.40 *T aw , and 0.20*T aw wall temperature conditions 
for slow wave acoustic disturbances are 9.47in, 9.67in, 14.62in, and 19.9in. respectively. Transition location is 
delayed with wall cooling for the same slow wave frequency because of the stabilization of first mode. However the 
amplitudes of wall pressure fluctuations increased with wall cooling to larger values than of adiabatic case. The 
delay of the transition locations and comparison of the amplitude levels can be easily seen in Figure 16.Figure 17 
shows the variation of eigenvalues ( r , (q) for the slow and fast waves for F= 1.2x15 4 at T w =0.20*T aw . 

Figure 18 compares the evolution of unsteady density fluctuations obtained from the simulations for slow wave 
(F=1 .2x1 5 4 ) at a fixed time for adiabatic and cooled wall (T w =0.20*T ai 7) conditions. Figure 18 (a) and (b) clearly 
show that there is no entropy layer generated for small bluntness case and the disturbances excite the boundary layer 
up to the wall. These figures clearly show the disturbance evolution in the nose region. Acoustic disturbances pass 
the bow shock and directly enter the boundary layer. Figure 16 (b) shows that the boundary layer is thinner and 
amplitude of the fluctuations inside the boundary layer is weaker for cooled wall compared to adiabatic wall 
presented in Figure 18(a). 

Figure 19 and 20 shows the propagation of density fluctuations inside boundary layer from leading edge to the 
end of the computational domain. In Figure 19(a) disturbances interact with the bow shock in the nose region and 
directly enter to the boundary layer. Disturbances form rope shape structure inside the boundary layer and on the 
oblique shock. The interaction is obvious between the density fluctuations on the oblique shock and boundary layer. 
Figure 19 (b) shows the disturbance field from 0.3in. to 2.0 in. In this figure four different zones are observed 
similar to previous studies 6, n " 12 . First zone is the area outside the shock where acoustic disturbances propagate 
uniformly. In the second zone acoustic waves transmitted through shock layer. The third zone is the area between 
the shock and the boundary layer. This region consists of transmitted external acoustic field and the disturbances 
radiated from the boundary layer. It is interesting to see that the third region is much quieter compared to the 
acoustic waves outside the shock layer. This implies that the acoustic waves are weakly transmitted through the 
shock. The fourth zone is the boundary layer where the boundary layer disturbances evolve. Figure 19(c) clearly 
shows decay of density fluctuation amplitudes from order of Iff 5 to 1 (I 8 which corresponds to the first mode 
stabilization. In Figure 20(a) and (b) density fluctuations inside the boundary layer are not discernable and the 
region between shock layer and cone wall seems quiet. However, acoustic disturbances propagated from 4 in. to 16 
in. in the order of 15 7 or below then started to gain amplitude and became visible (also look Fig. 15 (d)). This 
phenomenon shows the necessity of high resolution grid and higher order accurate schemes. Finally Figure 20 (c) 
shows destabilization of the second mode disturbances near the end of the computational domain. 

Discussion and Conclusion 

The receptivity and the stability of hypersonic boundary layers over a blunt cone with 5 degrees half-angle with 
nose radius 0.001 in. are numerically investigated at a free stream Mach number of 6.0 and at a Reynolds number of 
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7.8x1 (f /ft. Both steady and unsteady solutions are obtained by solving compressible Navier-Stokes equations using 
the 5 th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using a 3 rd 
order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. Unsteady flow is forced using 
slow and fast acoustic disturbance waves with the most amplified disturbances calculated based on mean flow 
analyses and also forced with non-dimensional frequency of F= 1.2x1 Of 4 . Computations are performed for different 
wall temperatures 804 (adiabatic wall), 603 ', 402, 322, 241 , and 161 °R. 

Table 4 summarizes computation parameters and the receptivity coefficients for the most amplified slow and fast 
mode acoustic disturbances at different wall temperature conditions. It is shown that the first mode of slow wave 
acoustic disturbances stabilized by wall cooling. However, wall cooling also caused the destabilization of the second 
mode and transition location moved to upstream from X=19.1 in for adiabatic case to X=7.61 in highly cooled wall 
case. It is also interesting to observe that wall cooling did not affect the first mode of fast acoustic disturbance 
waves. The receptivity coefficient of the fast wave case is 50 times greater than of slow wave case. We can conclude 
that the boundary layer is much more receptive to fast acoustic waves as compared to the slow wave. 


Table 4. Computation parameters and receptivity coefficients for the most amplified frequencies at different 
wall temperatures. 


F x 10" 


^ wall) 


wall' n V 




^ wall) 1 


0.96 s LOO 8.5x1 O' 6 4.23**’ # 194 0.552 

1.08 s 0.75 4.8 x 10" 5 2.39*’ # 13.1 0.0734 + 

1.32 s 0.50 6.2 x 10" 6 0.309* 10.7 0.0041 + 

2.0 s 0.20 2.0 x 10" 8 0.010** 7.61 0.0002 

2.0 s 0.20 I.OxlO -6 0.500** 7.59 0.0029 


T m = 113.98 °R, *P ac - 2 x 10 -5 , **P ac = 2 x 10 -6 . * These values divided by 10 to make the initial amplitudes 
in the same order. * in these simulation neutral point is not observed and amplitude values from x= 1 in. is used. 
s Slow wave, F Fast Wave. 

The transition locations stayed almost same for adiabatic wall and cooled walls (603 °R and 402° Ft) respectively 
for 9.475 in., 9.714 in. and 9.672 in. Transition locations increased dramatically for wall temperatures 322°R, 241 °R 
and 161 °R to 14.628, 17.013, and 19.906 inches. This is happened due to the stabilization of the first mode 
disturbances. However, amplitude of wall pressure fluctuations of cooled wall case ( T w = 161°R) increased 2.34 times 
of the adiabatic wall case ( T w = T aw =804°R). Table 5 summarizes the simulation parameters and gives the receptivity 
coefficients for the same slow wave disturbance frequency, F= 1.2x1 O' 4 , at different wall temperature conditions. 


Table 5. Computation parameters and receptivity coefficients for the same slow wave disturbance frequency 
at different wall temperatures. 



^ = 113.98 °R, P ac = 2 x 10 5 , Nose Radius R n = 0.001 in. # Neutral point is not observed for this case and 
amplitude from x= 0.207 is used. 

Because of the initial growth of first mode in adiabatic and cooled wall (T w =0.75*T aw ) neutral points are not 
observed. Therefore wall pressure fluctuations at 0.207 in. are used in the receptivity coefficient calculations. 
Neutral point locations moved downstream with wall cooling. Receptivity coefficients are 1.5225, 1.4613, 0.3246, 
0.0497, 0.0059, and 0.0023 respectively for wall temperatures 804, 603, 402, 322, 241, and 161 °R. 
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Wall cooling reduced receptivity coefficients and increased the transition Reynolds numbers. The receptivity 
coefficient for adiabatic wall case is 1.5225 and they are much smaller in the order of 10" 3 for highly cooled cones 
(241, and 161 °R). This raises some questions 34 about the transition process over cone with small bluntness. If the 
receptivity coefficients are very small for the second modes as was found in this paper, how can the amplitude of the 
disturbances reach such high values? Most probably non-zero acoustic incident angles may produce larger 
receptivity coefficients than zero incident angle used in this paper. 
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Figure 1. Sketch of the computational domain. 



and (c) Re-C r planes. 
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Figure 3. N-Factor curves for decreasing frequencies for a blunt cone (r n =0.001 in.) at different wall temperature 
conditions, (a) Adiabatic wall (T aw ), (b) T w =T aw *0.75, (c) T w =T aw *0.50, and (d) T w =T aw *0.20 
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Figure 4. Contours of temperature for flow over a cone with different surface temperatures at M=6.0 and Re=7.8x10 6 . (a) 
Adiabatic wall (T aw ), (b) T w =T aw *0.75, (c) T w =T aw *0.50, and (d) T w =T a w *0.20 




Figure 5. Density contours for flow over a cone with different Figure 6. Mach contours for flow over a cone with different 
surface temperatures at M=6.0 and Re=7.8x10 6 . (a) Adiabatic surface temperatures at M=6.0 and Re=7.8x10 6 . (a) Adiabatic 
wall (T aw ), (b) T w =T aw *0.20 wall (T aw ) f (b) T w =T aw *0.20 
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Figure 7. (a) Wall to free stream temperature ratio along the cone surface, (b) Effect of wall cooling on edge 
Mach number. 






Figure 8. Mean flow density profiles at different X locations in similarity coordinates for different wall 
temperatures, (a) Adiabatic wall (T aw ), (b) T w =T aw *0.50, (c) T w =T aw *0.30, and (d) T w =T aw *0.20. 
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Figure 9. Mean flow temperature profiles at different X locations in similarity coordinates for different wall 
temperatures, (a) Adiabatic wall (T aw ), (b) T w =T aw *0.75, (c) T w =T aw *0.50, (d) T w =T aw *0.40, (e) T w =T aw *0.30, 
and (f) T w =T aw *0.20. 
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Figure 10. Mean flow density profiles at different X 
locations in physical coordinates for (a) Adiabatic 
wall (T aw ), and (b) T w =T aw *0.20. 
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Figure 11. Mean flow temperature profiles at 
different X locations in physical coordinates for (a) 
Adiabatic wall (T aw ), and (b) T w =T aw *0.20. 


17 

American Institute of Aeronautics and Astronautics 


0.4 

- r 0 = 0.001 in. F=0.80*10' 4 Slow wave 

- 

0.2 

Pwa/P 00 

...... iiiiiiilll 

- 

0 

-0.2 



- 

-0.4 


- 


0 5 10 X (in.) 15 20 

(a) Adiabatic wall (T aw ) 


0.5 

i 0 0 

' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' r- 

Wall Pressure Fluctuations 

T wa „ = Cooled Wall (603 R) i 

F= 1 . 08 x 1 O ' 4 Slow Wave ,1 


; .iiiiiiilll 



Q. 

- 0.5 

: ^ 1(1 



1 : i i 1 : : : 1 : : : 1 : : : 1 : : : 1 : 



■' °0 2 4 6 8 10 12 14 


X (in.) 


(b) T w =T aw *0.75 





Figure 12. Pressure fluctuations on the wall, (a) Adiabatic 
wall (T aw ), (b) T w =T aw *0.75, (c) T w =T aw *0.50, (d) 
T w =T aw *0.20, and (e) T w =T aw *0.20 (Fast Wave). 
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Figure 13. Pressure fluctuations on the wall in Log Scale, 
(a) Adiabatic wall (T aw ), (b) T w =T aw *0.75, (c) T w =T aw *0.50, 
(d) T w =T aw *0.20, and (e) T w =T aw *0.20 (Fast Wave). 
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Figure 14. Pressure fluctuations on the wall under the 
effect of same forcing frequency F = 1 .2*1 0" 4 , (a) Adiabatic 
wall (T aw ), (b) T w =T aw *0.50, (c) T w =T aw *0.40, and (d) 
T w =T aw *0.20. 


Figure 15. Pressure fluctuations in log scale on the wall 
under the effect of same forcing frequency F = 1 .2*1 0" 4 , (a) 
Adiabatic wall (T aw ), (b) T w =T aw *0.50, (c) T w =T aw *0.40, 
and (d) T w =T aw *0.20. 
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Figure 16. Comparison of wall pressure fluctuations under the same forcing frequency (F = 1.2*10" 4 ) at 

different wall temperature conditions. 


F = 2*1 0 4 , T w = 0.20 * T aw = 1 61 R 



Figure 17. Eigenvalues (a r , on) for fast and slow modes for F = 2.0*1 O' 4 . 
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(b) Cooled wall (T w =T aw *0.20) 

Figure 18. Contours of unsteady density fluctuations inside the boundary layer near the nose region, 
(a) Adiabatic wall (T w =T aw ), (b) Cooled wall (T w =T aw *0.20) 
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Figure 19. Contours of unsteady density fluctuations inside the boundary layer along the cooled cone 
wall (T w =T aw *0.20), (a) x=[0.0-0.3]in, (a) x=[0.3-2.0]in, and (a) x=[2.0-3.7]in. 
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Figure 20. Contours of unsteady density fluctuations inside the boundary layer along the cooled cone 
wall (T w =T aw *0.20), (a) x=[4.0-5.7]in, (a) x=[6.0-15.0]in, and (a) x=[14.0-19.8]in. 
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